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We study Coulomb interactions in the finite bias response of Mach-Zehnder interferometers which 
exploit CO- propagating edge states in the integer quantum Hall effect. Here, interactions are partic- 
ularly important since the coherent coupling of edge channels is due to a resonant mechanism which 
is spoiled by inelastic processes. We find that interactions yield a saturation, as a function of bias 
voltage, of the period-averaged interferometer current which gives rise to unusual features, such as 
negative differential conductance, enhancement of the visibility of the current, non-bounded or even 
diverging visibility of the differential conductance. 
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Introduction. — Topological edge states in the integer 
quantum Hall effect (IQHE) [I] represent an ideal play- 
ground for testing coherence of electronic systems at a 
fundamental level. Harnessing of edge states as quasi 
one-dimensional (ID) chiral electronic waveguides has 
allowed the successful realization of a number of elec- 
tronic interferometric setups, such as Mach-Zehnder's 
[2H5] and Fabry- Perot-type's [6]. These devices ex- 
ploit counter-propagating edge states localized at op- 
posite sides of a Hall bar, which are brought in con- 
tact and mixed by quantum point constrictions (QPCs) 
that mimic the effect of optical Beam-Splitters (BSs). 
In Mach-Zehnder interferometer (MZI) setups the chiral- 
ity of electron propagation makes it necessary the adop- 
tion of non simply-connected, Corbino-like geometries, 
which limit the flexibility of the devices. In these sys- 
tems electron-electron (e-e) interactions are in general 
responsible for dephasing via the coupling with external 
edge channels [71 18], which manifests as a reduction of 
the visibility of the Aharonov-Bohm (AB) oscillations as 
a function of the bias voltage. In particular, puzzling 
behaviors have been reported [3H5] in the finite bias re- 
sponse of MZIs, in which the visibility presents a lobe-like 
structure and phase-rigidity [H [9Hl5] . 

An alternative MZI scheme inducing coherent mixing 
between edge states co-propagating at the same bound- 
ary of the sample, has been suggested [TBI as a more flex- 
ible architecture, which allows multiple device concate- 
nation [ini |T7j and represents an ideal candidate for im- 
plementation of dual-rail quantum computation schemes 
[l8l [To] , As schematized in Fig. [T] in such setup the 
BS transformations are now implemented through series 
of top gates, organized in arrays of periodicity A tuned 
to compensate for the difference between the momenta 
of the co-propagating channels (the inner channel i and 
the outer channel o of the figure). The corresponding 
AB phase difference is introduced instead by separat- 
ing the edge modes in the region between the two BSs 
through the action of a central, unmodulated, top gate 
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FIG. I; Sketch of the MZI. Two sets of A'^ top gates arranged 
in arrays with periodicity A and separated by the distance d 
represent the L and R BSs of the interferometer. A central top 
gate G locally lowers the filling factor to v = 1 and separates 
the two edge states (i and o), which experience a path length 
difference Ad and acquire an AB phase 4>ab. 



- see Refs. [111 [IHl IH] for details. The effects of e-e 
interactions in this MZI are likely to be qualitatively dif- 
ferent from those observed in the Corbino-like settings of 
Refs. [2HS]- Indeed, while in the latter the direct coupling 
between the channels that enter the interferometer can be 
neglected in the regions where they are coherently mixed 
(the QPCs), in the scheme of Fig. [l] this is no longer pos- 
sible due to the non negligible spatial extension of the 
top gate arrays. This implies a strong interplay between 
coherent mixing and interactions which might impair the 
MZI response. Aim of the present work is to target such 
interplay. As detailed in the following, we show that 
the inter-edge current / measured at the output of the 
setup of Fig. [l] possesses a strong non-linear dependence 
in the bias voltage V that, while still exhibiting AB oscil- 
lations, leads to saturation of the associated mean value 
averaged over the AB phases. As a consequence the de- 
vice presents negative differential conductance with un- 
bounded visibility. Such anomalous behavior occurs since 
interactions enable inelastic scattering which spoils, for 
increasing voltages, the coherence needed for inter-edge 
coupling to occur at the BSs. Furthermore, as long as the 
interactions between the edge channels can be neglected 
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in the region between the two BSs, we also observe that 
for large voltages the visibility of current gets enhanced 
with respect to the non-interacting case. 

Model: — Before discussing the role of e-e interac- 
tions in the device of Fig. [l] we find it useful to briefly 
review the basic properties of the scheme in the non- 
interacting case. The underlying idea |20) is to imple- 
ment BS transformations between two co-propagating 
channels {i and o), via the action of a pair of arrays of 
top gates (see Fig. [I]) which are spatially modulated at 
periodicity A. Following Refs. [20l [21] we describe them 
through potentials of the form t[^{x) = sin^(7ra;/A), 
for — AA^ < X < 0, and tji{x) ~ tflsin^(7r(a; — d)/X), 
for d < X < d + \N {N being the number of elements 
of a single array while d being the distance between the 
two sets as measured according to the coordinates of the 
channel o). Introducing the difference Afc = ki — ko be- 
tween the momenta ki and k^ of the two edges, the tun- 
neling amplitude at a given BS can then be expressed 
as ic<sinc(AA^(Afc — 27r/A)/2), at lowest orders in (here 
a = L,R indicates the left and right BS, respectively) [20]. 
In this scenario Mach-Zhender interferences can be ob- 
served by introducing between the two BSs a top gate 
which locally lowers the filling factor to = 1 and di- 
verts the inner edge state toward the interior of the mesa 
[IB] . This way, the two channels are guided along paths of 
difference lengths, do — d and di — d+Ad, thus acquiring 
an AB phase difference 0ab proportional to the magnetic 
field B and to the area enclosed by the path and a dynam- 
ical phase. The transmission probability T(e) at energy e 
of the MZI, from inner channel on the left to outer chan- 
nel on the right, is then given by T{e) oc [I^lP + I^^Rp + 
2|iLti?| cos(0AB + eAd/ui?)]S'^, where vp is the group ve- 
locity of edge states and S = sinc(AA^(Afc — 27r/A)/2) 
is the BS amplitude which is optimal when the resonant 
condition Ak — 2tt/X is met. A bias voltage V applied 
between channels i and o gives rise to a zero-temperature 
current I(V) oc J^^ de T(e), whose associated visibility 
V/ = (max^J — min0j)/(max0/ -|- min^/) of the AB os- 
cillations amounts to V/ = VCT*'''|sinc(eV^A(i/2i'F)|, with 
Vi°^ = 2|iL<fl;|/(|tL|2-h|ifl|2) being the oscillation visibil- 
ity of the differential conductance a = dl/dV. 

To analyze the effect of interactions in such a set up 
we describe the linearly dispersing electronic excitations 
around the Fermi energy by means of two chiral fermion 
fields e**'''"^Vm(a;), with m = i,o, each propagating at 
mean momentum km- The kinetic Hamiltonian can then 
be written as {h = 1) TJkin = -™Fj2mI dxipl^d^ipm- 
The action of the BSs are instead assigned by means of 
the tunneling Hamiltonian iJtun = J2a=L hi^a + ^a), 
where Aa describe the action of the gate arrays potentials 
and are defined as Al = J dx tL{x)e^^'^^ ^pl{x)^pi{x) , 
Ar = Jdx tfl(a;)e'^'=^ iljl{x)ip,{x + Ad) with ii(x), 
tfi(x), and Ad introduced previously. In these expres- 
sions the local phase shift e*^*"'^ accounts for the reso- 



nant behavior of the MZI. The AB phase of the setup 
(together with a dynamical phase term kidi — kgdo) is 
instead included in the tunneling amplitude of the right 
beam-splitter, i.e. tn — >■ e^'^^^tji. 

As far as e-e interactions are concerned, an electron 
propagating in one edge channel can interact with all the 
electrons in the Fermi seas of both channels. Although 
the precise form of the interaction potential is unknown, 
screening provided by top gates makes it is reasonable 
to assume a short range density-density interaction. The 
latter however needs not to be uniform in the whole sam- 
ple: as a matter of fact, while intra-channel couplings 
are present everywhere, the inter-channel interactions de- 
pend on the edge channel spatial separation, which in our 
setup varies strongly from place to place (see Fig. [T]) . In 
particular in the region between the BSs edge states are 
brought far apart by the central top gate G and it is rea- 
sonable to assume the inter-channel coupling to be off. 
Inter-channel interactions, however, cannot be excluded 
in the regions before and after G, where the tunneling 
term i?tun is active. Indicating with p„i = '0m'0m the 
ID density operator in channel m = z, o, we account for 
these effects by introducing an inter-channel e-e coupling 
^^intcr = S dx J dx' : po{x)U {x , x') pi{x') : characterized 
by a coordinate dependent potential U{x, x') which nul- 
lifies in the central top gate region (i.e. < x < d) and 
which approaches the short range behaviors 2'iTg5{x — x') 
and 27TgS(x — x' — Ad) on the Ihs part (i.e. x < 0) 
and on the rhs part [x > d) of the setup, respectively 
(the transition between these regions being smooth). In 
these expressions the interaction strength g has the di- 
mension of a velocity while, similarly to fftun, the pa- 
rameter Ad accounts for the relative coordinate shift ex- 
perienced by the inner edge with respect to the outer. 
A short range intra-channel coupling term of the form 
-ffintra = '^^J2m I ^x Pmi^) ■ ^^^o Considered (in this 
case however no coordinate dependence is assumed). 

Methods: — Including the e-e interaction, the response 
of the interferometer driven out from equilibrium by a 
bias voltage pi — fi^ = applied between the two edge 
states, will be analyzed at lowest order in the tunnel- 
ing Hamiltonian i/tun (an exact analytical treatment be- 
ing impossible). At second order in the amplitudes ta 
this allows to express the current flowing through the 
outer edge as I{V) = e^^^^ / di([Aj^(t), A^(0)]) , where 
Aa{t) — 6^^°*^ Aae~^^°* is the tunneling term evolved 
through the kinetic and interaction components of the 
system Hamiltonian, i.e. Hq = i/kin + -^intra + -Winter, 
while expectation values are taken with respect to the 
ground state of the system biased by the chemical poten- 
tial difference eV [7]l22]. 

Expanding the summation over a and /3, we recognize 
the presence of three contributions: I — II + Ir + 
with /a = e / dt{[A'l{t),Ac,{0)]), for a = L,R being di- 
rect terms which do not depend upon the relative phase 
accumulated by the electrons when traveling through the 
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MZI, and with /$ = e J dt{[A\{t), Ar{Q)]) + c.c. being a 
cross-term, which is sensitive to the AB phase. Exphcit 
expressions for these quantities are obtained by means 
of the two-point electron and hole Green's functions 
gUx,t;x') = {^mix,t)i;lix',0)) and g'^ix,t;x') = 
{tp}^{x, t)tlj„i{x' , 0)), which we calculated by bosonization 
of the Hamiltonian Hq. For this purpose we introduce 
chiral bosonic fields 0„i, which satisfy the commutation 
rules [<^m(a;), 0m' (a;')] — iTrS„i^m'Sign{x — x') [ISllllj and 
express the fermion fields as ibm = ^ii=e^^*^'"^/^e~*'^'" 

v27ra 

where a is a cutoff length that regularizes the theory at 
short wavelengths, are the Klein operators, — 
J dx : Pm{x) '■ are the total number operators of the edge 
states, while finally L is the edge quantization lengths. 
Observing that pm = (l/27r)9x0m + Nm/L, the kinetic 
and intra-channel Hamiltonian (apart for an irrelevant 
term proportional to J2m ^rn) becomes -ffkin 

the renormalized edge group velocity [71 and 
where 7?™*''^ — T^vJ2,n^m/L is a capacitive contribu- 
tion. Vice-versa, exploiting the smooth variation as- 
sumption of the potential U{x,x') and the fact that L 
is the largest length in the problem, the inter-channel 
interaction term yields 
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{d.Ao{x)){dxd^^{x + Ad)) + H'S'"', (1) 



with i/^''^'^ — {2Trg/L)NoNi being a cross capacitive con- 
tribution. The Hamiltonian Hq can now be brought 
to a diagonal form by solving the eigenvector equation 
[^^o,7±(e)] + e 7±(e) = 0, which defines bosonic energy 
eigenmodes 7±(e). Accordingly we obtain 



4>mix,t) = 



de 



s=± 



f'mi.x.,e)ls{e)+^-C. 



(2) 

where the wavefuctions ip'^ satisfy the relation 

Es/o°° f Im[¥'m(a^,e)¥'m'(a;',e)1 = f (5™,™'Sign(x - x'), 
and 1/r > is an energy cutoff related to a via the 
non-interacting dispersion e — vpq, i.e. a/r = vp. 
Solving the equations of motion by requiring only con- 
tinuity of the wave functions we find ip^^ = f±{x) for 
X < 0, (p^ ^ C±e*'^^/'' for < a; < d™, with C+ = 1, 



C_ = 0, and (^^ 



3"'*"-/V±(a; - d™) for x > d™. 



in terms of the symmetric and antisymmetric combi- 
nations f±{x) = (e*^^/^+ ± e*'^/'"-)/2. Two new ve- 
locities enter the problem, a fast charged mode which 
propagates at u+ = v + g and a slow neutral mode 
which propagates at V- = v — g. Exploiting these re- 
sults, the electron and hole Green's functions can be fi- 
nally written as = e-*'^"(*-(^-"^')/''-'^;„(a;, t; x') and 
gh^ = e^^'"(*-("^-^')/"-)^;™(a:,t;a;'), with the zero-bias 
Green's function Qm{x,t; x') = 



1 /gi0™(£c,t)g-i0„(a;',O)\ 



which, thanks to the quadratic nature of the bosonized 
Hamiltonian Hq, can be easily computed in terms of the 
wave functions tpf^. In particular to compute the direct 
current terms 1^ and Ipi we only need the correlation 
functions in the BS regions: x,x' < for the left BS, 
and X, x' > dm for the right BS. At zero temperature, for 
these combinations we find 
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(3) 



(here t = t—ir). The term is obtained instead through 
crossed combinations. In particular for x > dm and x' < 
we get 
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2vrt;F Vi^ - ^')/vs + dmil/v - - t 

(4) 

and analogously for x' > dm and a: < 0, with the replace- 
ment dm ^dm- 

Currents: — With the help of Eqs. ^ and Q the 
I-V characteristics of the setup can now be explicitly 
computed. In particular for the direct contributions of 
the current one gets 
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(5) 



with the resonance function S{x) = sinc(A7V(AA; — 
27r/A)/2 — a::/2), and the density of states at the Fermi en- 
ergy np = XN/{2vp). The energy scale Cc = {XN{l/v- — 
l/v+))~^ is associated to the difference in time-of-flight 
for propagation of the bosonic excitations at speeds v+ 
and V- through the BS of length AiV: in the absence of 
interactions it diverges. We notice that while for a low 
bias la is linear in V, for a bias larger than ec it satu- 
rates to the constant value 1^^^ = 7rGon'^\ia\'^ec. This 
is shown by the dashed lines in Fig. [2] top panel, where 
/„ is plotted as a function of V for different values of g. 
Such behavior is due to the fact that inelastic processes, 
which are induced by the interaction and increase with 
increasing voltage bias, spoil the resonant coherent effect 
which is responsible for the transfer of charges between 
the two edge channels in the BS, thus suppressing Jq,. 

The cross term contribution /$ to the current is 
obtained instead through integration over the branch 
cuts ^2j of the Green's functions Q. Introducing 
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(6) 

together with the phase (py(e) ~ e/€c + {Ak — 2TT/X)XN + 
eV[Ad/{2v^) + {d + Ad/2){l/v - l/v_)], we obtain 

I^iy) = ^n%2\iLtR\\A\ cos(0AB + arg(^)), (7) 

where Jo[x] is the Bessel function of the first kind and 
Edyn = V / Ad. The latter is a new energy scale associated 
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to the dynamical phase difference acquired by the elec- 
trons in the interference region between the BSs, where 
inter-channel interaction is absent and excitations move 
at speed v. Due to the bias dependence of arg(y^), the 
current /$ shows oscillations versus bias around zero. 

The overall I-V characteristics (full lines plotted for 
different values of g in Fig. [2j top panel) show an 
oscillating behavior which becomes non-monotonic for 
large values of the bias V. The associated visibility is 
Vi{V) = V<°^|yt|//;^de S'2(e/eJ, where vi°\ as in 
the absence of interactions, provides the upper bound 
(V/ < Vi°^ < 1). Plots of Vi{V) for different values of g 
are shown in the bottom panel of Fig. [2] For g/vp = 0.1 
(red curve), weakly interacting regime, Vi{V) closely fol- 
lows the non-interacting case (black curve) for V <V* — 
2TThvp/{eAd) and thereafter oscillating without reaching 
zero. In the strongly interacting regime g ^ vp, the scale 
€c already dominates at relatively low biases, leading to 
a completely different behavior. As shown by the blue 
curve {g/vp = 10), Vi{V) decreases with the bias very 
rapidly up to eV < tc and very slowly thereafter (with 
the parameters used in the plot ec/{&V*) ~ 0.32). The 
almost constant value exhibited by the curve is the con- 
sequence of the fact that the current oscillates around 
a constant value, with amplitude of oscillations that de- 
creases very slowly with bias. Interestingly, in this regime 
the visibility in the interacting case is higher than in the 
non-interacting case. Very peculiar is also the case of a 
symmetric interferometer. Ad — 0, where edyn diverges 
and the visibility becomes V^''^^ = vi°V2(^)l/-^o(V"), 
where Fp{V) = j^^''''^^''^ dx sinc^{x). We find that 

the ratio Vf^'*"°/vi-°^ goes, for eV ^ Cc, asymptotically 
to the finite value 21n(2)/7r. Therefore, a symmetric de- 
vice over-performs with respect to an asymmetric one for 
large voltages. 

The highly non-monotonic behavior of I{V) detailed 
above is best understood by studying the differential con- 
ductance dl/dV. In particular, its direct contribution 



amounts to 



dV 



°^np\ia\ S {eV/tc), and it shows how 
e-e interactions effectively give rise to an interaction- and 
bias-dependent shift of the resonance condition 
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As a result, dIa/dV becomes negligible beyond an ap- 
plied voltage on order of Ec, which explains the saturation 
of the period-averaged current in Fig. [2] in the strongly 
interacting regime (notice also that for weak interaction 
the applied bias voltage V allows for fine tuning of the 
resonance condition). At the same time a non- negligible 
cross term dl^/dV results in an overall behavior of the 
total differential conductance dl/dV which exhibits re- 
gions of negative values. Furthermore, while the associ- 
ated visibility of the AB oscillations in the presence of in- 
teractions is known to exhibit values larger than one |12j , 




FIG. 2: (Color online) Response of the MZI versus bias V 
[in units of V* = 2nhvF / (e Ad)] for different values of the 
interaction parameter g/vp- We have assumed u = g, 0AB = 
0, I^lI = ~ \t\ and that the resonant condition, Afc = 
2tt/\, is satisfied. Top panel: total current / (full curves) and 
BS current (dashed curves) in units of /* = 2nj,\t\'^e^V/h. 
Bottom panel: Visibility of the AB oscillations in the current 
Vi in units of Vi°\ According to Ref. [20], we set XN/Ad = 3 
and d/Ad = 1. 



for our setup this quantity diverges for values of bias V 
such that the direct contributions to the differential con- 
ductance vanishes. 

Conclusions. — The analysis presented in this paper 
confirms for a MZI with co-propagating edge states the 
picture in which electron-electron interaction in general 
reduces the performances of an electronic interferometer. 
There are, however, striking differences with respect to 
non-simply connected MZI architectures. In the latter 
case interactions couple the interferometer edge channels 
to additional channels that carry information away from 
the system [SHSl [12] . In our setup interactions directly 
couple the two interferometer edge states and informa- 
tion is redistributed in the system. The major impact of 
interaction effects is the spoiling of the resonant single- 
particle tunneling, that realizes the BS, while the oscil- 
lating component of the interferometer current is only 
marginally affected because inter-edge interactions are 
negligible in the interference region. This leads to the 
unexpected result that strong interactions yield a reduc- 
tion of the current visibility for small voltages, but an 
enhancement for larger voltages, with respect to the non- 
interacting case. Furthermore, the differential conduc- 
tance can become negative in some voltage range, while 
its visibility can take large values or even diverge. 
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